Quantum Stochastic Processes: A Case Study 



Michel Bauer *'* and Denis Bernard * 

* Institut de Physique Theorique de Saclay^, 
CEA-Saclay, 91191 Gif-sur-Yvette, France. 

* Laboratoire de Physique Theorique de I'Ecole Normale Superieure, 
CNRS/ENS, Ecole Normale Superieure, 24 rue Lhomond, 75005 Paris, France 



Abstract 

We present a detailed study of a simple quantum stochastic process, 
the quantum phase space Brownian motion, which we obtain as the 
Markovian limit of a simple model of open quantum system. We show 
that this physical description of the process allows us to specify and 
to construct the dilation of the quantum dynamical maps, including 
conditional quantum expectations. The quantum phase space Brow- 
nian motion possesses many properties similar to that of the classical 
Brownian motion, notably its increments are independent and iden- 
tically distributed. Possible applications to dissipative phenomena in 
the quantum Hall effect are suggested. 
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1 Introduction and motivations 



Quantum Brownian motion [T] is usually defined by coupling quantum me- 
chanically a massive particle to a bath of harmonic oscillators, see e.g. 
refs.[2l[3lllj and references therein. This is not the quantum analogue of what 
probabilists call Brownian motion (which is the continuum limit of random 
walks) but that of a Newtonian particle submitted to friction and random 
force (which is also called Brownian motion in the physics literature) , whose 
dynamics is described by mx + ax = f with a the friction coefficient and / 
the random forcing. These two notions are related: the former is the over- 
damped limit, or the massless limit, of the latter. The mass plays a role of 
short distance cut-off, and the velocity x scales as 1/ ^/rn. It diverges in the 
probabilistic Brownian motion which is known to be nowhere differentiable. 
This causes troubles when defining the quantum analogue of the probabilis- 
tic Brownian motion as the massless limit of the quantum Brownian motion 
because the momentum p = mx scales as p ~ y^, whose massless limit is 
incompatible with the canonical commutation relation [x,p] = ih. The first 
aim of the present paper is to understand how to consistently formulate the 
over-damped quantum Brownian motion. 

A possible way to overcome these difficulties may consist in going one 
dimension higher and consider the two dimensional Brownian motion. An 
advantage is that one may then couple the Newtonian particle to an out- 
of-plane external magnetic field B, so that the classical equation of motion 
becomes 

mz + (a + ieB)z = f 

in complex coordinate z = x + iy, with e the electric charge. If the force 
/ is white-noice in time, i.e. {f{t)f{s)) oc 6{t — s), the massless limit still 
describes the probabilistic 2D Brownian motion, and its quantification will 
provide a candidat for the 2D (over-damped) quantum Brownian motion. 
The magnetic field only rotates and dilates the Brownian trajectories. Since 
these trajectories are known to be statistically conformally invariant |6j, the 
effect of the magnetic field is classically quite innocent. Quantum mechan- 
ically the situation is different. In absence of random forcing, the energy 
spectrum is that of the Landau levels. The magnetic length is = y^h/eB 
provides a short-distance cut-off which would allow us to take the massless 
limit. Since the cyclotron frequency ojb = eB/m diverges as m vanishes 
(or as B increases) , taking the massless (or the over-damped) limit amounts 
to project on the first Landau level. Thus a model candidate for the over- 
damped 2D quantum Brownian motion is a particle confined to the first 
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Landau level and coupled to a harmonic oscillator bath|^ Recall that, once 
projected onto the first Landau level, the two coordinate operators x and y 
are not commuting [9| but form a pair of canonical operators, [x, y] = so 
that the motion actually takes place in the quantum phase plane. The model 
obtained by following this projection strategy leads to a coupling between 
the oscillator bath and the particle which (although quadratic) is yet too 
complicated Hoping, or claiming, for universality property, we consider 
an alternative model with a much simpler coupling between the quantum 
phase plane and a bath of harmonic oscillators, see eq.Q. 

In this way, in a weak-coupling/long-time/continuous-density limit, we 
obtain a Markovian quantum process, which we call "quantum phase Brow- 
nian motion" and which bares many similarities with classical Brownian mo- 
tion. A similar, but not identical, model has been considered by E. Davies 
in ref.[5] but without analysing the quantum process. 

The two dimensional Brownian motion possesses very peculiar geomet- 
rical properties linked to conformal invariance, see e.g. ref.[7j. It would be 
very interesting to know how these properties are deformed after quantisa- 
tion. To understand the quantum geometry of the quantum 2D Brownian 
motion is another long term motivation of the present work [8] . 

The second aim of the paper is to provide a description ~ as complete as 
possible - of the quantum stochastic process defined by this simple model. 
We aim at presenting the algebraic structure underlying the process which 
goes beyond the dynamical maps and the associated master equation, which 
code for the evolution of the reduced density matrix |im [TH [T^ . This is in 
particular needed if one is willing to evaluate, or even define, multi-point 
expectations. Recall that a classical stochastic process is defined over a 
probability space equipped with a probability measure and a filtration (i.e. 
an increasing family of cr-algebras) to which the process is adapted. The 
process is specified by a consistent set of finite-dimensional distributions 
(encoded in transition kernels for Markov processes) and the filtration al- 
lows to define conditional expectations. The filtration encodes the increase 
of knowledge as time and the process go on. In Appendix [6] we recall the 
algebraic structure induced by these data. Quantum mechanically, transi- 
tion kernels are replaced by completely positive dynamical maps acting on 
the algebra of observables of the quantum system (but not of the bath). By 
duality, this specifies the evolution of the system density matrix and it is 

■^One may wonder whether this model may have relevance to the description of dissi- 
pative effects in the integer quantum Hall effect, although it does not deal with the edge 
currents. 

^We nevertheless hope to report on this model in the near future [5|. 
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enough for evaluating one-point functions (but not multi-point). Most, but 
not all, of the physics literature only deals with density matrices. This is 
however not enough to fully specify the process and in particular to define 
the quantum analogue of conditional expectations. This requires extending 
the algebra of observables and the dynamical maps to a set of embedded al- 
gebras (quantum analogues of the filtration, which roughly keeps track of the 
process up to time t) with flows intertwining them and projectors (quantum 
analogues of conditional expectations) with compatibility properties. This 
extended algebraic structure, which mimics that of the classical processes 
(see Appendix [6]) , is called a dilation of the dynamical maps [lOl [TTl [12l [13] . 
Given completely positive maps, there is generically not a unique dilation 
but a set of compatible dilations. We present this structure in the simple 
(but quite generic) model we describe and show that the physical description 
of the process given by coupling the quantum system to an oscillator bath 
specifies the dilation]^ 

Exemple of dilations of dynamical maps have of course been given in 
the literature, see e.g. refs. jlOl [TH [T2] and references therein. Elements 
of the construction we present here also bare some similarities with the 
quantum stochastic calculus of ref.[16]. However, as far as we know, our 
construction of the dilation from the microscopic physical model has not 
been described earlier. Our model provides a simple framework in which 
concepts of quantum stochastic processes may be explained and exemplified. 

The paper is organised as follows: In the section [2| we define the model 
and describe the long-time/weak-coupling limit in which it becomes Marko- 
vian. In section |3] we describe the quantum process with a detailed presen- 
tation of the quantum filtration, of the dynamical maps and the associated 
Lindbladian, of the flows and the conditional expectations, etc... The quan- 
tum process we have defined enjoys many properties analogous to those of 
Brownian motion. In particular we show that its increments are independent 
and identically distributed and we show that the conditional expectations 
are neutral with respect to both left and right mesurable multiplication^ 
We also make contact with processes introduced earlier on the mathemat- 
ical literature. In Appendix [6] we recall the algebraic structure underlying 
classical stochastic processes and in Appendix [7] we present an alternative 
approach to the weak-coupling/long-time limit. 



^We show this property in the case of this simple quadratic model, but we hope to also 
proof it in a future publication |8] in the case of an arbitrary finite dimensional quantum 
system linearly coupled to a bath. 

^This is explained precisely below. 



5 



2 A simple heat bath model 



2.1 The model 

Out-of-equilibrium quantum systems may often be described by coupling 
them to reservoirs which, quantum mechanically, may be modelled as col- 
lections of harmonic oscillators. Here we consider one of the simplest model: 
a harmonic oscillator (which we call the system oscillator) linearly coupled 
to a reservoir set of independent harmonic oscillators (which we call the 
bath). We denote by e the energy of the system oscillator and by Ua the 
pulsations of the reservoir oscillators. The system creation-annihilation op- 
erators will be denoted w'^ and w and those of the reservoir by and a^. 
They satisfy canonical commutation relations: 

The other commutators vanish: [aa,aa'\ = and [i(;^,aa] = ['ii),aa] = 0. 

The hamiltonian of the system oscillator plus the reservoir is the sum 
of three pieces: H = Hb + )^Hsb + Hs where Hb is the bath hamiltonian, 
Hs the system oscillator hamiltonian and Hsb a coupling hamiltonian. Ex- 
plicitely 

We set /i = 1 and A is a dimensionless controlling parameter and are 
coupling constants. The Hilbert space is the tensor product of the Fock 
spaces associated to each of the pairs of creation-annihilation operators: 

H = rs^rB with Fb = (g)„ r„. 

At initial time, the system is supposed to be prepared with a density 
matrix pQ^ps where the system density matrix po is yet unspecified. Density 
matrices are normalised, Trp = 1. The bath density matrix is choosen to 
be gaussian, pB = 0„ Pa with pa = e'^" and = 1 — e~^«. We 
shall parametrize the bath density matrix via the mean occupation numbers 
Ua = Trr^ipaolaa) with = l/(e'^" - 1). 

Expectations of product of observables (which may be evaluated at dif- 
ferent times) arc by definition given by their traces weighted by the density 
matrix: (O) = Tr-^(po "X" Pb O). As usual with composite systems, we first 
trace over the reservoir Hilbert space as we are not aiming at observing it. 
We denote by E the trace over the reservoir: 

nO] = TvTs{pBO). 



6 



By abuse of language we shall call E the expectation, it maps general oper- 
ators into operators acting on the system Hilbert space r^. The (complete) 
expectations are thus given by (O) = Trr^ipo^iO]). 

2.2 Weak-coupling/long-time behaviour 

The aim of this section is to describe the model in the weak-coupling/long- 
time/continous-density limit, corresponding to A — )• with A^t fixed, in 
which the dynamics is known to become Markovian \14:\ IT5] . The continuous- 
density limit ensures dissipation because the bath relaxation time scales 
as the inverse of the bath energy level spacing. Computations presented 
in this Section are a bit technical but they are needed to get an explicit 
description of this limiting behaviour. The output formulas are summarised 
. They are going to be the starting point of the construction of the 
quantum stochastic process which we call quantum phase space Brownian 
motion. Computations may be skipped (in a first reading) if the reader is 
only interested in the output. 

In the Heisenberg picture, the equation of motion for an operator O is: 

O = i[H, O]. 

In our model, this leads to linear equations for the oscillators: 

w + iew 

da + iuJa0.a 

(together with their adjoints). Our aim is to understand these equations 
in the weak coupling limit for an idealized bath involving a continuum of 
frequencies. This can be done routinely using Laplace transform. 

Define the Laplace transform / of a function / by / = /q^°^ /(t)e~P*. 
Then the equations of motion become: 

pw{p) - w{0) + iew{p) = -iXy^ipadg^p) 

a 

pda{p) - aa{0) + iuJadaip) = -iXCpaWip). 

We solve the second equation for da and reinject in the first to get 



m 



eqs.(3p 



-iX}^(paC 

a 

-iX (paW 
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As we shall discuss in details in Appendix^ the function p+ze+A^ \^£^lL_ 



dp ^p, ^(o)-^AE.^«»(o) 

2m p + ie + A2 



can only vanish for purely imaginary values of p. So Laplace inversion tells 
that 

rc—ioo 

where c is an arbitrary positive constant. 

In the weak coupling limit the exchanges between the bath and the 
system take a long time, and the short time evolution of the system is 
dominated by the energy scale e. Once this is factored out, the relevant times 
are longer. As we shall see below, if the bath is idealized by a continuum of 
oscillators, the relevant long time scale will be of order l/A^ and there is a 
non-trivial zero-coupling limit for 

z{t) = w{t/X^)e'''/^\ 

Let us define 

^ X^p -ie + lOJa 



After translating and rescaling p in the above formula for w{t), we get a 
formula for z{t) which is a sum of a term proportional to the initial value 
2;(0) and terms proportional to aa(0), that is: 

z{t) = f{t)z{Q) + with m = E /„(t)a„(0). 



where 



C—lOO 



dp ^pt 



2m p + r{p) ' 



' C— JOO 

J c—ioo 



2i-K [X'^p — ie + iuJa] [p + r(p)] 
Taking complex conjugates, we note that 

Jc-ioo 227r [X^q + ie- lUJa] [q + T{q)\ 

Because the bath is Gaussian and ^{t) is linear in bath oscillators, its 
properties are fully specified by the value of the commutator 
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and the two point function E[^^(s)^(t)]. Prom the basic oscillator commu- 
tation relations, wc have — /a(*)/a(^)- Using the identity 



_2 i_ _ 

p+u q-U} p+q 



'Vq{^ + ^) we get 



J c—ic 



dp dq ^pt^g^ _ r(p) + T{q) 
2m 2m \p + T {p)\[q + T {q)\[p + q] 



In the same vein, we find E[^^'(s)^(t)] = '^^riafaifyfais), which can be 
written as: 



2m 2m \p + T{p)] [q + T{q) []p + q] 

where 

Turning now to the situation when there is a density of oscillators, we 
replace |(^o!p<^(^ ~ ^a) by r(a;) and, interpreting as n(ua), we find: 

+°° r{u) _ /■+°° , n{uj)r{uj) 



X^p — ie + 



Taking naively the limit — >• 0"*" and remembering that the integration 
contours stay in the half- plane Jiep > 0, we find that r(p) 7 with 

7 = f(0+) = 7rr(e) + i/^°°da;^^. (2) 
Jo e — uj 

Analogously, T(p) + T{q) — > 27rn(e)r(e). 
Thus, in the weak coupling limit we get: 



"°° dp ^ 1 



2m P + 7 ' 



rc—too 



Note that E[^t(5)^(i)] = ri{e) [m,^Hs)]- 
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The remaining integrals are obtained by application of the residue theo- 
rem. For instance, f{t) = e~'^^. The commutator requires a slight discussion. 
If t > s > 0, 



-ps 



-■ys 



^ 2i7r q + jp + q 1-p 

because the integration contour can be pushed to the left. This has no pole 
at p = 7. Then, 



-ps 



2iTT p + 7 1 — P 



-■ys 



-7(i-s) 



-ys—yt 



7 + 7 



because the integration contour can be pushed to the left again. The case 
s > t > is treated analogously. 



Hence, defining 
Git,s) 



^i{t—s)^m y 



we end up with f{t) = e and 

z{t) 



g-|i-s|5Re7 _ ^-{t+s)?Re-y 



Git,s) 
n{e)G{t,s). 



The above derivation of the long-time/weak-coupling/continuous-density 
limit relies heavily on the magic of contour deformation. Though the final 
result, damping when the bath is infinite, is physically satisfactory, it may 
seem disturbing that taking the limit of a function with many poles but all 
purely imaginary, one finds functions with single poles having a negative real 
part. The model we consider is simple enough that one can make a more 
down-to-earth derivation which shows more clearly that damping arises from 
destructive interferences. It also leads to an explicit spectral representation 



(21 ) of the process directly inherited from the discrete harmonic bath oscil- 



lator decomposition. This is presented in Appendix [7[ 



3 Quantum phase space Brownian motion 

Let us reformulate (and summarise) the output of the weak-coupling/long- 
time limit of previous section. The time evolutions of the system canonical 
operators z and z^ are: 



z{t) = e-^'z + m , z^ (t) = e-"z^ + e (t) 



-It J 



(3) 



10 



with 7 a complex parameter, 5Re 7 > 0. We shall set 7 = K+iv and 7 = K—iv. 
The fields ^(t) and which represent "the quantum noice", are linear 

combinations of the reservoir creation- annihilation operators. Commutation 
relations are: 

[z,zt] = i 

and 

The system operators commute with the reservoir operators so that [z, = 

The measure E, induced by tracing over the reservoir degrees of freedom 
with the bath density matrix, is such that i\t) and ^(s) are gaussian with 
two-point function: 

m\t)i{s)\=rioG{s,t) 



with no real positive. In previous Section \2?2\ this parameter was denoted 
n(e). Since physical phenomena concentrate at the energy scale e in the 
long-time/ weak-coupling limit, all details of the reservoir parameters have 
been erased and summarised in the very few parameters 7, no and in the 
form of the commutation relations G{t, s) and in the two-point functions 
noG(s,t). 

The aim of this section is to extract (and illustrate) the structure of 
quantum processes in this particular example which possesses all generic 
properties plus some peculiar ones. It has most properties for being a quan- 
tum analogue of a (two-dimensional) Brownian motion. In particular its 
increments are independent and identically distributed. 



3.1 Flow and expectations 

• Algebras: We first have to identify where the flow takes place. Let ^0 be 
the algebra generated by the system canonical operators z and z\ and let 
i?[o,t] be the algebra generated by all the ^(s) and i\s) for s G [0, t]. We set 

At = Aq® B^Q^t] , Aoo = Ao®S[o,oo)- (5) 

Although understandable by "common sense" and useful in practice, the 
definition of i?[o,t] "tbe algebra generated by the C(s)'s with s < t" is 
not precise. Thus, let L2(M+) be the Hilbert space of functions s — )• f{s), 
such that e'^^ f{s) is square integrable on the positive real line, equipped 
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with the scalar product {g\f)^ = ds g{s) 2Ke'^'^' f{s). Let r(L2(M+)) 
be the Fock space over L^(M+), and let QJj be the canonical operators on 
r(L^(M+)) depending linearly on functions / with commutation relations 
[23/,5jJ] = (9|/).. Then 

More explicitely, the bath operators ^ (t) may be represented as (recall that 
7 = K + iu): 

^{t) = e-^' [ du{2Ke^^'')^ a{u), 
Jo 

where a{s) and a^(s) are bare canonical operators with commutation rela- 
tions: 

[a(s), a{s')] = 0, [a{s), a^(s')] = <5(s — s'). 

One may use this representation of ^(t) to write a "quantum stochastic 
differential equation" for z{t). Indeed differentiating eq.Q, we get: 

dz{t) = -^z{t)dt + V2Ke'''^d%{t), (6) 

where we set d^{t) = j/^'^* a{u)du. A similar representation is of course at 
the basis of the quantum stochastic calculus due to Hudson and Parthasarathy 

m- 

Let Boo = -S(r(L^(M+)) be the algebra of bounded operators on r(L^(M+), 
which we identify with B[o,oo)- Recall that r(L2(M+) = r(L2([0,t]) ® 
r(L^(]t, c«)). We may then define -B[o,t] as the i?oo-subalgebra (this is more 
an embedding than an inclusion): 

5[o,t]=i?(r(L2([0,t]))cSoo. 

Notice that -B[o,s] ^ for s < t and B^o = -B[o,t] (8> B^f^^o)- There is some 
freedom in the choice of the algebras B^q^^, the above choice is a minimal 
one, maybe not adapted to all situations one may encounter. However, most 
of the following computations are algebraic and valid for a larger class of 
operators on r(L^(M4.)) than those in -B[o,t]- 

• Flow: The time evolution of the system operator defines a flow Jt from 
Aq to At by: 

Jt-.Ao^At, Jt{z) = zit) = e-^'z + m- (7) 
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Similarly, Jt(z^) = z^{t) = e + ■^^(t). It is a *-liomomorphism since 
[z{t),z^t)] = 1, because G{t,t) = 1 - e-(T+^)*, so that 

with 

Eq.Q is the quantum analogue of a stochastic classical flow with ^(t) playing 
the role of quantum noise. 

• Measure: The quantum noises ^{t) are elements of B^o- The measure over 
it is E. It is gaussian with two-point function E[^^(t)^(s)] = noG(,s,t) so 
that 

^[^^^l£.Ht,)^^^,atl) . . .e/^iv«t(t^)gMivC{t^)] = exp[^ /ii//jG(ti, tj)(no + l{i<i})]- 

Here /ij and /ij are formal parameters not necessary complex conjugate. If 
one chooses = —jJ-p then the operators involved in the above equation 
are in Boo- 

The measure E is neutral with respect to elements of the system oscillator 
algebra, that is E[ao b] = aoE[6] for any € Aq, so that it extends to a map 
from Aoo to Aq: 

E : Aoo ^ Aq. 

On the bare canonical operators a{s) the measure is E[a^(s)a(s')] = no(5(s — 
s'), and: 

E[d^{t) d2l(t)] = 0, E[d2lt(t) d2t(t)] = no dt. 

To complete the setting, one has to provide the measure Trr^ {po • • • ) 
on Aq specified by the system density matrix. Classically, E would be the 
measure on the noise and Trpg {po • • • ) the measure on the initial position of 
the stochastic process. 

3.2 Dynamical maps and density matrices 

• Dynamical map: The map <I>t on is obtained by evaluating one-point 
functions: <I»t(a) = E[Jt(a)] for a G Aq. In the present case, it yields: 

Of course $j(l) = 1. As can be checked by a direct computation, it defines 
a semi- group on Aq: 

$t o $5 = 
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The form of the generator of the dynamical maps that we are going to find, 
in eq.Q below, implies that are completely positive maps, as required 
by the theory of quantum processes [TOl |TT1 |T2] . It is the quantum analogue 
of transition probability kernels for classical stochastic processes. 

• Lindbladian: The Lindbladian is the generator of the dynamical semi- 
group. Indeed, the semi-group law (with extra continuous properties) im- 
plies that 

$t = exp( tL ) 

where L = ^^t\t=o is a so-called Lindbladian, with a structure imposed by 
the general theory of completely positive dynamical maps |17| : 

L{a) = i[hs, a] + XI C^DjaDj - DjD^a - aDjD^.), a e Aq 
J 

with some effective hamiltonian hg = hi and some operators Dj and gj > 0. 
In the present case, L is quadratic in z and and, for all a ^ Aq, 

L{a) = —iv \z^ z, a] + g\ {2zaz^ — zz^a — azz^) + §2 (2z^o2 — z^za — az^ z), (9) 

where we set 7 = «; + zz^ and 

gi = Krio, 52 = K(no + 1). 

Of course L(l) = and L{a'') = L(a)"l'. Set L{a) = —ii'[z'^z, a] +D{a) where 
D{a) is the so-called dissipative part of the Linbladian. It may alternatively 
be written as 

D{a) = — Kno {\z\ [z, a]] + [z, [z\ a]]) + k {2z'^ az — z'^ za — az^ z). 

The double commutator term coincides with the (quantum analogue of the) 
two dimensional Laplacian. As expected, the imaginary part of 7 only enters 
into the hamiltonian part of the Linbladian, but note that the frequency u 
is not the bare frequency e. The dissipative part of L only depends on Ke 7 
and riQ. This Lindbladian has been considered in ref.[3] to describe damped 
harmonic oscillator. 

• Master equation: This refers to the evolution equation for the system 
reduced density matrix. It is obtained from the observable time evolution 
by duality since Trr3(po ^t{o)) = Trrg(/9t a) for a £ Aq. Hence pt = e*^* -po- 
So that 

^Pt = L*.pt (10) 
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with 



L* ■ p = iu[z^ z, p] + §2 {2zpz^ — pz'^ z — z^ zp) + g\ (2z^ pz — pzz^ — zz'^ p). 



Eq.(lO) is the quantum analogue of Fokker-Planck equations for Markov 
processes. It preserves the normalisation of p, i.e. lh:Yg{L* • p) = 0. 

The system admits an invariant measure which is a density matrix pinv 
such that L* • piaw = 0. It is gaussian: 

Pin. = \e-''''\ withe'^ = '^^, 
z no 

and Z = 1/(1 — e~"). Alternatively, no = 1/(6°^ — 1). An effective tempera- 
ture may be defined either by cj = u/T or a = e/T. The invariant measure 
depends on no but not on k. 

Since k represents damping effects, the approach to equilibrium is k- 
dependent. It is exponentially fast as e"^*^* with relaxation time 1/k. It is 
governed by the smallest eigenvalue of L* which is —2k: 

L* ■ pi = -2k pi, 

with eigenvector pi = e~'^^^^[z^z — 2no]. Other eigenvalues of L* (in the zero 
charge sector) are —2pK with p integer and eigenvector pp = e~'^ ^'^ ^ Pp{z^ z) 
with Pp polynomial of degree p. 

The invariant measure is stable since all eigenvalues of L* are non posi- 
tive. This is a consequence of the complete positivity of the dynamical maps 
which implies that the dissipative part oi L*-p is the sum of terms of the form 
gj{2DjpD'^-— pD'^-Dj—D'^-Djp) which are all non positive operators for gj > 0. 

Indeed Trr^ [p {L* •/?)) is the sum of 2gj (Trr^ {pDjpD^)- Trp^ {pDjD^p)) < 
by the Cauchy-Schwarz inequality. 

3.3 Quantum filtration and conditional expectations 

• Filtration: The set of embedded algebras At = Aq® B[o,t] defines a filtra- 
tion 

Aq C As C At C Aoo, for < s < t < oo. 

The algebra As has to be understood as the past of the process up to time 
s. The image Js{Aq) of the system algebra by the flow is identified as the 
present of the process at time s: Js{Aq) C As. With these definitions the 



*We do not enter here in the subtleties associated to left/right time limits as the 
construction is clearly time continuous. 
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past includes the present. The algebra (g) -Bj^^oo) may be identified with 
the futur of the process. 

• Conditional expectations: These are projectors from into Ag such 
that 

Es-.Aoo^As, E^-^ o = Ejnin(3^_52), for Si, §2 >0, (11) 

and compatible with the flow and the stochastic map in the sense that, for 
all a € Aq: 

Es[Jt{a)] = Js[<^t-s{a)] , fort>s>0. (12) 
Es has to be understood as the expectation conditioned on the past of the 



process up to time s. Eq.(12) expresses the Markov property of the process. 



In the present case, compatibility with the flow imposes that, for t > s, 

]E_,[gM?''"(t)gA?(t) ] = gM(t-s)C'l"(s)g;i(t-s)^(s) ^noltip- fJ.it- s) p.{t- s)] ^ 

This defines E<j on single-time operator but we have to define it on product 
of multi-time operators. Since we know the commutation relations in A^q 
it is enough to define it on product of time ordered operators. That is, we 
have to define 

]g^[gMi5'''{*i)gA«i?{*i) . . . gA'JV?'''(tiv)g/i]V?(i]v) j 

with ti < ■ ■ ■ < ti\i. A constructive way to do it consists in (recursively) 
imposing that E^ is neutral with respect to left multiplication by elements 
of As, i.e. Es[a6] = aEs[6] for a G Ag. After a few computations summarised 
in Appendix [8j we get: 

]g^[gA«i5'''{*i)gAi?(*i) . . . gMJvC''"(t]v)gA«]v€(*Jv) j (^i^^ 

for s < ti < • • • < t]y, with 

^s^^ = ^fiUJ-jih-A) -^Mti;s)lJ'jitj;s) 

i<j i<j 



where tj-i = tj — ti and tj-s = tj — s. Notice that the r.h.s. of eq.(13) only 
involves ^(s) and ■^^(s). One checks that Es=o = E as it should be. 
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• Markov and other properties:. The previous definition is designed to ensure 
that relations (11) and (12) hold, but this can be checked directly. See 
Appendix [8} 

Since it only involves ^(s) and C^(s), it also satisfies a (weak) Markov 
property because the image by Kg of the future algebra -B]s,oo) maps into the 
present algebra Js{Aq): 

Es : 5],,^) ^ JsiAo) (14) 

By construction, is neutral with respect to left multiplication by ele- 
ments in As, but remarkably it is also neutral with respect to right multi- 
plication. Namely, 

Es[ab] = aEs[b] , Es[ba] = Es[b] a , for a G A, (15) 

for all b € Aoo ■ The first relation is true by construction and the second one 
follows from the relation 

^^mo),cHt)] = ^^it-s) mo),^\s)], 

for to < s < t. See AppendixjSj This property is remarkable. In particular it 
means that we would have got the same conditional expectations if we would 
have constructed it using neutrality under right multiplication instead of 
neutrality under left multiplication as we did. In this sense, the conditional 
expectation is unique. 

By construction (also), the conditional expectation satisfies the nested 
formula, similar to that valid in the classical theory, but for time ordered 
products s < ti < ■ ■ ■ < t]\f: 

Es[Jtj{ai) ■ ■ ■ Jtj^_^iaN-i)JtNiaN)] (16) 

for ttj £ Aq. Note that the nesting is ordered from right to left. A similar 
nested formula applies for an opposite ordering of the time but with an 
opposite nesting. Both formula follow from eqs.(ll) and from the left and 
right neutralities. 

The nested formula eq.(16) is of course an echo that the bath dynamics 
has been frozen by the infinite- volume/continuous-density limit. 

• Martingales: As for classical processes, we may define martingales ~ which, 
for instance, are instrumental for computing probabilities of events related 
to stopping times. A family of elements Mt in At is called martingale if 

Es[Mt]=Ms, fors<t. 
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Examples are for instance given by: 

To first order in p,, this gives e'^*^(t) as a martingale. 
3.4 Time evolution in Aoo 

• Flow: By eq.Q we have defined the flow of system observables. We may 
extend this flow to elements of ^oo- That is, we may define a semi-group of 
^oo-endomorphisms, as with fio = Id, such that: 

as - Aoo ^ Aoo, crsoat = as+t, fx^ o = J^+g, (17) 

for s,t > 0. The second relation, namely as o Jt = Jt+s, means that as 
extends the flow J^. In particular as acts like Js on elements of Aq. It 
implies that: 

a.(z) = e-^^z + e(s), as{m)^at + s)-e-'^'as)- (18) 

Similarly as{z^) = Js{z^) and as{i^{t)) = ^t(t + _ e-^*£}{s). One may 
check by a direct computation that as are *-endomorphisms and 

Ws{m).^s{£}{m = [m.i\t')] = G{t,t'), 

for all t, t' > 0. Of course, as o at = as+t- Since both as and Jt are 
morphisms, the compatibility relation as o Jt = Jt+s is valid on Aq because 
it is true for the generators of Aq. 

• Covariance: The flow at is compatible with the conditional expectations 
in the sense that 

atoEs = Et+s o at. (19) 

This is a direct consequence of the previous construction. Indeed, neutrality 
of for elements of As implies that this is true on Aoo if it is true when 
applied to elements of the future algebra -B]s,oo)- To check it on B^g^oo), 
it is enough to check that it holds true for product elements of the form 
Jti{o-i) ■ ■ ■ Jt^iaN) with ttj S Aq and times ordered s < ti < ■ ■ ■ < tjy. On 



such elements, eq.(19) is a direct consequence of the nested formula eq.(16) 
for conditional expectations and of the compatibility relation at{Js{cL)) = 
Jt+sia) on ylo- 
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3.5 Independence of increments 

We call (Ts(^(t)) = S,{t + s) — e~^*^(s) the increment from s to t + s. 

Since ag is an endomorphism and since the two-point functions are pro- 
portional to the commutators G{t,t'), the increments are identically dis- 
tributed: 

E[asi^Ht))aMt'))] = mHt)m]=noG{t',t) 
for all s > 0. 

Furthermore, the relation [^{to),^^{t + s)] = fi{t) [^{to), ^H^)] to < s 
and t > 0, which was instrumental in proving the left/right neutrality of the 
conditional expectations, implies that: 

Mm) , ^\to)] = 0, forto<s (20) 
n<^s{m)^Hto)] = 0, forto<s 

These two relations then imply that [o"s(^(t)), a] = and E[(Ts(^(t)) a] = 
for any a £ Ag. It means that the increments commute with the past algebra 
As and are independent of the past. Another remarkable property. 

4 Limiting cases and decoherence 
4.1 The quantum t-Brownian limit 

The process simplifies in the limit 7 — )• 0. Set 'j = k + iiy and let = 
^^(f)e~''^*. In the limit k — )• 0, these can be written as xi^) = fodsa{s), 

in terms of the bare canonical operators a(s), with [a{s), a^{s')] = d{s — s'). 
Then 

[x{t),xHs)] = 2min(t,s) 
nx\t)x{s)] = 2nomm{t,s) 

in the limit k — )• 0. The real and imaginary parts of xit) firs two Brownian 
motions but which do not commute except at time t = 0. 

This coincides with the quantum process (also called a quantum Brown- 
ian motion) defined in ref.[18] using Hopf algebra techniques. Actually, our 
initial algebra Ao has to be slightly modified in order to take care of the limit 
K — )• 0. Following ref.p2], we define the initial algebra as the Heisenberg 
algebra ^0 generated by elements z, and r with commutation relations 

\l,l'\=r, [r,zt] =0 = [r,5]. 
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The flow is defined by Jt{z) = x(^) and Jt(T) = t. It is a *-homomorphism. 
The dynamical maps is the 7—7-0 hmit of eq.Q, namely: 

In this limit, the conditional expectations are 

]E_,[g»?X'''Wg'?x{*)] — ^VxKs) ^rix{s) g2r?f?no(i-s)^ 

It indicates that the increments + ■5) ~ x{^) &re independent and identi- 
cally distributed. 

Note that the above limit is not the same as the limit k — )• but with 
D = Ktxo fixed which would gives 

[m,i\s)] = 0, m\ms)] = 2I?e-(*-^)min(t,s) 

and conditional expectations 

for t > s. It describes two commuting Brownian motions. The system 
algebra is still that of the harmonic oscillator, [z,z^ = 1. The fiow simply 
corresponds to add the random C-number ^(t) to the canonical operators 
z, an operation which obviously preserves the commutation relations. The 
Linbladian also simplifies in this limit, since only the double commutator 
remains. 

4.2 Semi-classical limit 

To take the semi-classical limit we put back the h factors and set e = hujQ 
and replace (pa by hipa so that fa has the dimension of a frequency. The 
hamiltonian is then: 

a a a 

The large time limit is taken in the same way. The factors h compensate in 
evaluating z{t) so that we again have z{t) = e'^h + Cit) with [C{t),CHs)] = 
G{t, s) where 

G{t, S) = e-^*-T^ (g{7+7)min{t,s) _ 

The two point function is still given by E[^"''(t)^(s)] = noG{s,t). The 
occupation number no may be estimated as if it is given by thermal equilib- 
rium at some temperature T, that is no = l/{e^^^/'^ — 1) and no — T/Hujq 
as h goes to 0. 
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The semi-classical limit is thus ^ — t- with noh = T/ujq fixed. 

Set 7 = K+iu and let X{t) = Vh^e z{t)e-'''^ and Y{t) = ^/^9m z{t)e-'''K 
Then, in the semi-classical limit, X{t) and Y{t) are two commuting gaussian 
processes with 

E[X{t)X{s)] =E[Y{t)Y{s)] = ^e-'^(t+«)(e2«min(M) _iV 

They are two commuting Ornstein-Uhlenbeck processes. As k — )• 0, X{t) / ^/k 
and Y(t)l ^[k goes over two (un-normalized) commuting Brownian motions. 



4.3 Decoherence and thermalisation 

Let us first consider the evolution of a density matrix which initially is simply 
the pure state /?o = with |q;) the coherent state \a) = e~°"/^ e"^^ |0). 

At time t, the density matrix is determined by duality via Trpg {pt a) = 
Trrs(po^i(a)) with given in eq.Q. Using Trrs(/Oo e^'^^e'^'^) = e°'^+'^^, 
this gives 

Tirgipt e^^'' e^^) = e"^(*)+"'^(*) ^no{^l^^-^J,{t)^^{t)) ^ 
This is solved in terms of superposition of coherent states by: 

Pt= \a) e "t"- ^ " {a\. 

J T^fTt 

with at = no(l — e^^*^*) and a{t) = e~^*a, an evolution equation similar 
to that of z but without the quantum noice. The interpretation is clear: 
the density matrix is an (incoherent) sum of coherent states centred around 
the un- noisy solution a{t) with dispersion at as in the classical case. As it 
should be, it converges at large time towards the thermal invariant measure, 
since: 

f dada _m 1 -az^^z 

Poo = / \a) e "0 a = — e ^ 

J Trno no + 1 

with 6°^ = (no + l)/no. The thermalisation time is 

''"therm = 

Let us now consider an initial density matrix which is still that of a 
pure state but for a state linear combination of coherent states, a kind of 
Schrodinger cat: 

po = {u\a) + v\^)){u{a\+v{^\), 



®With our convention, X has dimension of \/h since X ~ v^5Re z. To go to the 
physical picture with X having dimension of a length, one would have to introduce the 
mass scale m of the harmonic oscillator and set X ~ -^/fi/mcjo 5fte z. 
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with |a— /3| large. We explicitely write po = uu Pq°'+uv Pq^+vu Pq'^+vv p^^ . 
By linearity of the evolution equation, the density matrix at time t will be 
of form 

- aa I - a/3 . - Ba , - SB 
pt = UU Pf + uv p^ + vu p^ + vv p^ . 

The evolution of each of the terms is determined as above. For instance: 

ft = (a|p) / \a) e "t {a\, 

J vat 

or equivalently, 

(6|pf"|c) = {a\P){h\c) e-^('^-^W)(^-^W), 

with 1 6) and |c) coherent states. In order to code for the decoherence, let 
us consider the ratio of the matrix elements of the 'mixed' density matrices, 
pf" or p"'^, over the pure ones, p"" or p^'^, defined by: 

^ (6|pf |c) {h\pf\c) 



{h\pr\c) (6|pf|c)' 

with 9^0 = 1- An explicit computation shows that this ratio is independent 
of 6 and c and given by 

Recall that |(a|/3)p = e-l"-'^l^ At short time, at ~ 2no«;t, so that 

log^Ht ~ -2(no + l)Kt\a - /5|^ 

at short time. Hence, the mixed density matrices decrease exponentially 
faster than the pure ones. 

Thus there are two regimes: first a short decoherence time Tcoher = 
K(no+i)|a-/9p ' whose semi- classical limit is (recall that tiq 2± T/hujQ and 5X ~ 

^coher 



kT \5X\^' 

and which decreases as the inverse power of the square distance between the 
states, and then a longer thermalisation time rtherm — l/n- 
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5 Conclusion 



Motivated by an initial problem related to the quantum Hall effect, these 
notes have investigated a simple open quantum system. Starting from a 
traditional coupled system-reservoir model, we were able to study most of its 
salient features in the weak-coupling/long-time limit, from the more physical 
to the more mathematical. 

The model is simple enough to be solved by several methods. One of 
these concentrates on the limiting behavior of expectations and commuta- 
tors, without paying attention to the limit of operators themselves. The 
output is an abstract (almost axiomatic) representation result for the lim- 
iting process in Fock spaces. Another approach, which leads to the same 
result, pays more attention to the relationship between the bath oscillators 
of the initial model and their weak-coupling/long-time cousins. As expected, 
dissipation arises via destructive interferences. 

Along the way, we retrieve a number of well-studied objects that we 
connect together, like the damped harmonic oscillator on the physical side, 
or the quantum stochastic calculus and quantum geometry on the mathe- 
matical side. We pay a particular attention to the explicit construction, in 
this simple case, of a quantum filtration and the associated quantum con- 
ditional expectations. This leads to the specification of the (mathematical) 
dilation of the quantum dynamical maps. It gives access to multi-time cor- 
relation functions that cannot be computed via the knowledge of the system 
time dependent density matrix only. The quantum phase space Brownian 
motion we obtain possesses peculiar properties analogous to those of the 
classical one. We observe the remarkable fact that despite non-trivial non- 
commutative effects, the conditional expectations have the property to be 
transparent to multiplication by measurable observables (on the left or on 
the right). We also do a sample computation of decoherence. 

One of the missing probabilistic concepts is that of stopping time - which 
is instrumental for computing exit probabilities or harmonic measures. The 
concept is defined in the mathematical literature, but the examples are few 
and a bit artificial. In the present case, and if the system starts in the 
ground state, one can define the time at which the system leaves the ground 
state. This can be checked to be a stopping time, i.e. deciding if it occurs 
before time t can be done by looking only at conditional expectations up 
to time t. Alas, due to continuous time, this time vanishes almost surely. 
However, our discussion of filtration and conditional expectations can be 
extended straightforwardly to a discrete version of eq.Q, namely zq = z 
and Zn = Xzn-i + (1 — lAp)-*^/^ a„ for n > 1 (where |A| G [0, 1[ and the a^s 
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are standard independent annihilation operators). Then one checks without 
difficulty that, starting the system in the ground state at n = 0, the time at 
which it leaves the ground state is a stopping time, which follows a geometric 
distribution. Other examples of meaningful stopping times, in the discrete 
and especially in the continuous time case, would surely be desirable. 

We hope to come back in a near future to all these themes in other sim- 
ple open quantum systems - including generalisations to finite dimensional 
quantum systems and the associated quantum stochastic processes, aspects 
of quantum geometry of quantum phase space Brownian motion - and to the 
more puzzling case of dissipative effects in the quantum Hall effect and their 
relations to the behaviour of the first Landau level coupled to a reservoir. 
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6 Appendix: Algebraisation of stochastic processes 



The aim of this appendix is to recap the algebraic structure underlying classi- 
cal stochastic processes for comparison with the quantum theory. Although 
we did not find it gathered in one reference, we have no doubt that the 
material presented here is already somewhere (dispersed) in the literature. 

Suppose we are given an adapted stochastic process, defined on a prob- 
ability space S7, with measure E and filtration Tt-, and valued in some mea- 
sured space V with measure dy. For a; G let t — >■ ^^(a;), or simply Xj, 
be the realisation of the process. By construction this defines a path on V 
(for simplicity we may assume that it is continuous). The measure E then 
induces a measure on paths on V . 

To simplify the writing, we shall restrict ourselves to time-homogeneous 
Markov processes but most of the algebraic structure we shall describe re- 
mains valid for non Markovian processes (up to simple reformulations). Let 
Pt{dy,x) be the transition kernel, i.e. the probability density for the proc- 
cess started at x at time t = to be in the neighbourhood of y at time t, so 
that for any (good enough) function / from V to, say, C we have: 

E,[/(Xt)] = J f{y)Pt{dy,x). 

Let ^0 be the algebra of measurable functions on V with values, say, in 
C (with the Borel a-algebra) , At be the algebra of Jt-measurable functions 
and some inductive limit of these algebras adapted to the filtration: 
Ao = Punc(y C) and 

At = {f& Func(n — >■ C), s.t. / is Jj-measurable}. 

Since to any a; G we associate a path t — )■ Xt{uj) on V, function of At are 
functions on parameterized paths [0,t] — >■ F of "length" t drawn on V. The 
algebra Aoo is the algebra of E-measurable functions. We identify Ao with 
At=o- Of course, since J^t is a filtration, we have the embeddings: 

AoCAsCAtC Aoo, for s <t. 

Naively, J^s-measurable functions are functions only of the portion of the 
path t — > Xt{uj) for t G [0, s], i.e. they code about the knowledge of the 
process from time to time s. 

The data of the stochastic flow Xt (tautologically) induces a map Jt from 
Ao to At which to any test function f on V associate an J^-measurable 
function by: 

Jt: A -> At, Jt{f){oj) = f{Xt{u)). 
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Given Pf, we may define the stochastic map on by: 

$i(/)(x) = E,[/(Xt)], feAo. 

Alternatively, $f = Eo J^, with a slight abuse of notation. As Pt, they define 
a semi-group ° = ^t+s- Note that $t is defined using only one-point 
functions (involving only one time). 

Conditional expectations are maps from At onto .4s for s < t: 

E[Ft\Ts]eAs, for Ft G A- 
They are neutral with respect to multiplication by J-'^-measurable functions: 

E[FsF\J^s] = FsE[F\J^s]. 

As a consequence, they are nested projectors since E[F5|Js] = Fg for Fg G .4s 
and 

E[E[F\Ts] l-^t] =E[F|J-^i,(,.,)], 

for any E-measurable function F. 

The data of Pt{dy, x) specifies the process. In particular, by conditioning 
recursively on the position of the path at intermediate times, multipoint 
correlation functions may be evaluated as 

E,[hiXt,)---fniXtJ] 

= J fiiyi)PtMyi,x)--- J fniyn)Ptn-tn-i{dyn,yn-i), 

for ti < ■ ■ ■ < tn- By definition conditioning on Xt=Q = x amounts to put 
the Dirac point measure on J^q, so that these multipoint expectations may 
also be written in terms of the maps $t as: 

EiMXt,) ■ ■ ■ fn{XtJ\To ] = ^t, (/l ^t,-H {■■■ fn-1 ^t„-*„_i(/n))) 

More intrinsically (and for non Markovian processes), using the nested 
property of conditional expectations, for any J^^^ -measurable functions Fj, 
we can recursively write 

E[Fi---F„|Jo] = E[Fi---F„_iE[F„|j;„_J|Jo] 

= E[ FiE[F2 • • • E[F„| J • • • I -Ft^ll-^o ] 

We see that this formula can be written in terms of maps $t;s, which map 
Fj-measurable functions onto J^-measurable functions, defined for s < t by: 

$t;s(F) = E[F\J^s] e As, for F G At- 
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Indeed, 



E[ Fi • • • F„| J-Q ] = cD.^.o (Fi ^trM {■ ' ' Fn-l (^n))) 

For Fj = fj[Xt-) and time homogeneous Markov processes, this last formula 
coincides with the previous one since then (^t\s{f{Xs)) = ^t-s{f) for s < t. 



7 Appendix: Another approach to the long time 
Umit 

The aim of this Appendix is to re-derive the long-time/weak-coupling/continuous- 



density limit of Section 2.2 differently in way which makes more transparent 
the fact that damping arises from of an infinite number of destructive inter- 
ferences. 

Without loss of generality, we assume that the ipaS are non vanishing. 
We also assume that the Wq-'s are all distinct. 
Define a rational function 

|2 



s{n) = n- e - x'^Yl 



a 



Q — OOn 



Its degree is equal to the number of oscillators in the bath plus one. Direct 
substitution shows that Q is an eigen-frequency of the coupled system (i.e. 
the coupled equations have a solution for which both w and the Oq's have a 
time dependence e~*^*) if and only if S{^1) = 0. The variations of S show 
immediately that the solutions are real and intertwine the w^'s: as expected, 
there is one more solution than the number of oscillators in the bath. In 
particular, there is at most one negative eigenfrequency. Let {^^k} be the 
set of eigenfrequencies of the coupled system. 

The Hamiltonian is bounded below if and only no solution is negative, 
i.e. if 

We assume that this condition is satisfied in what follows. 
Writing 

1 _ \ ^ Ftk 

s{n) ~ ^n-Qk' 
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Figure 1: The typical shape of the function S{Q). 



we claim that the general solution of the coupled system is 



k \ a 



a 



— I ^k — ^(3 



These formute could be obtained mechanically via residue calculus using 
the inverse Laplace transform solution presented in Section[2.2[ But a direct 



argument is easy too. For fixed k, the coefficient for 0/3(0) is just ^^^^^ times 
the coefficient for w, so this is indeed a family of solutions of the coupled 
equations. What remains to be checked is that the initial conditions are 
satisfied. This is the case because 



k 



k 



Yl 7?s ^ = '\2^ — ™y ^■ 



k 
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The first equation expresses that S{Q) behaves like Q for large Jl, so that 
1/S'(0) behaves Uke The second equation expresses that S{n) has a 

pole at each Ua so that 1/S{uja) = 0. For the third identity, we consider 

two cases. If a ^ B, writing 77^ k^, r = , ^ , ( o ^ , „ ^ . I , the 

third identity is a consequence of the second. For a = f3, the left-hand side 
is nothing but —^5^ evaluated at f2 = Wq. As S{Q) ^ ~^ J^" ^ for J7 
close to the result follows. 

From now on, our aim is to compute the behavior of the O^'s and R^s 
when the w^'s and ip^s are adjusted to mimic a smooth density. We choose 
a smooth positive density function r{uj) and consider Riemman sums with 
a mesh rj, a small parameter with the dimension of a frequency, leading to 



ze[i,L] 



n-rii' 



Later, we can let 77 — >■ 0"'" and L +00. We assume that e > X];e[i l] 
which is ensured for r] small enough if e > /q^°^ r{uj)^, an integral that we 
assume to be convergent. The zeroes of 5 can be ordered as < < ?7 < 
Hi < 2?7 < • • • < ^L-i < Lr) < Ul- Though not stricly evenly spaced, the 
fifc's are good enough to define Riemman sums. 

The phases of the (faS can be reabsorbed in the bath oscillators, so we 
may assume that the (paS are real and positive. Then, we can write the 
coupled time evolution as 



wit) = y: i?.e-^^^M-(o)+ E -m^^^] 

fce[o,L] V ie[i,L] ' ) 

Fix < -qL and take k such that krj <Vt < {k + l)r]. We want to find 
an approximate formula for 5(0) for fixed CI when 77 is small and L large. 
We shall see that the condition 77 log L small is important. 

More precisely, we let 77 but write = r]{k + 6), 6 g]0, 1[, so that the 
interval between two adjacent eigenvalues is blown up to size 1. We need to 

r{r]l) 
Aell,L] k-l+5 



compute r]{k + 5) — e — \^ Z^/efi l] k-l+s — ^- ^P^i* ^ 



rjrjl) ^ rjrjk) ^ / r(770 r{r]l) \ ^ rjrjl) 
^ k-l + S S ^ \k-l + S k-lj ^ k-V 
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The last sum is well approximated by 

- + 00 



Jo 



^ — OJ 



up to 0{r)). In the remaining sum, writing r{rjl) = r{rjk) + {r(r]l) — r{r]k)), 
the second term is 0(r? log L). When these errors can be neglected, a good 
approximation for S{fl) is 

\ tanvrc) Jq il — oj J 

This equation gives an approximation for the il-kS close to 0, by choosing 
5 = 6k such that the right-hand side vanishes. This is not directly needed, 
but it is the crucial ingredient to control in this region. Indeed, to get 
S'{Q) we can take the derivative with respect to 5 and divide by ry. This 
leads to r]S'{n) = X'^Tr^r{n) (l + ^^^). But S'{nk) = 1/Rk, leading to 

= ^ ^ - r^RiSi) 



for kij close to Q. So indeed, in this approximation is obtained by 
discretizing a smooth function. For instance. 



f+OO 

V RkC-'^"^^ / dne-'^^Rip). 
77^ Jo 



fce[o,L] 

This is enough to treat the term proportional to w{0) in w{t). To deal with 
the oscillator contribution in w{t), we need to control the sum 

k£[0,L] " ' 

By contruction, at t = this sum vanishes, and we may subtract = 
e-^^'^Efcelo.Ll^fen^ to get 

1 r+oo —iQt p—ir/lt 

fee[o,L] K I ju I 

If we define a/ = ^a{r]l) then in the continuum limit we get [a{uj),a^ {co')] = 
d{uj — uj'), and we find 

r+oo r+oo „—iQt „—iijJt 

w{t) = w dne-'^^R{n) + X dUdoj a{u) R{n) . 

Jo Jo il — LO 
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The limit A ^ 0+ in z{t) = w{t/\^)e^^*f^^ is now straigthforward. We 

set — 7- e + A^Jl, — 7- e + X^uj. Observe that Aa(e + A^w) is still a nor- 
malized oscillator, which we denote by a^^. Noting that R{e + A^fi) 
x'i(n^r{e)^+ln-f{e))^) ' ^hich we denote by Rq, we find that 

/+0O r+oo —iClt p—itot 
dne-'^^Rn+ / dndcj a^y^Ru 
oo J-oo il — u; 



i 



oo 



7r2r(e)2 + (17 - f(e))2 



-oo 



n-uj ''7r2r(e)2 + (J]-f(e))2- 



Integration over 12 is responsible for destructive interferences. Perform- 
ing this integral and using that 7 = 7rr(e) + if{e) leads to 

This gives the spectral representation of z{t) = ze~^* + ^{t). One checks 
that = G{t,s) as should be. 

The spectral representation also allows to check easily that the coupled 
bath-reservoir system is still hamiltonian, with Hamiltonian 

H = J dujual^auj + ^ ^ j duj (^z^a^ + za]^ + '^m^z'^ z 

Checking that this leads to the spectral representation is done routinely via 
the Laplace transform, interpreting at some point J duo as / du = 
TT for 3?ep > 0. Note that the hamiltonian is not bounded below due to 
the appearance of negative frequencies. This is simply due to the fact that 
taking the long time limit leads to shift the origin of energies by e. 



The relation between the coupling constant \J^^^ and the friction co- 
efficient 5fte 7 is no surprise (though it is usually interpreted as the Einstein 
relation only after introduction of a temperature) . 

To diagonalize the dynamics, define (recall that Rq = r(e)/|r2 -|- Z7p) 

\R^J V^ + n 7-0O + 

One checks that the 6's are normalized oscillators 
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such that [-ff, 6q] = —0,bfi, and that 

/ + 00 
-oo 

Just as z = dil R^'^bfi commutes with the operators and a^, one 
checks that duj a^^ is a normahzed anihilation operator that com- 

mutes with the operators 6q and 6q. 



8 Appendix: Proofs 



Here we present the computation which leads to the formula eq.(13) for the 



conditional expectations. Eq.(12) defines on single-time operator. We 



have to define it on product of time ordered operators, that is: 
]g^[gMi5'''(*i)gAi5{*i) . . . f^^^Ni\tN)Q^^Ni{tN) j 

with ti < ■ ■ ■ < t^. We impose that E^ is neutral with respect to left 
multiplication by elements of Ag, i.e. Es[a6] = aEs[6] for a £ Ag. As a 
consequence, it is enough to consider times bigger than s, that is s < ti < 
■ ■ ■ < t]y. We may then recursively compute the conditional expectation by 
imposing the required condition E^^ o E^j = Ejiiin(sj g^)- Indeed, inserting 
first Eg = Eg o ^tpf-i ^ < tN-i in the above conditional expectation and 
using that all operators e^^^^^^^^e'^^^^*-'') with j < N — 1 are left neutral with 
respect ^t^-i (by construction), we obtain 

]g^[gMiC'''(ti)g/iiC(ti) . . . gA«iv5U*iv)gMivC{tiv) j 

= Es[e'^^^^*^*^-'e'^^^^*^-' • • • e'^^-i^^(*^-i)e'^^"^^*^*'^-i^ E^^ [gMJv?''"(*]v)g/i]v€(tjv)j j 
— ]E^[e'^i?^(*i)gAi?(ti) . . . gA'iv-iC'''(<jv-i)gAjv-i?(t]v-i) 

^ gA']v(tjV;]V-l)C^(tjV-l)gMJv(*iV;JV-l)?(tjV-l) j gno[A«]V/ijV— MJv(t]V;JV-l)Ajv(*iV;JV-l)] 

with tj\f;N-i = tN — tN-1 and where in the last equality we used the formula 
for conditional expectations of single-time operators. The key point is that 
now the conditional expectation E^ only involves the — 1 times ti < • • • < 
tN-i, SO that we can recursively apply this procedure to fully compute the 
conditional expectation. After appropriate reordering of all terms, we get 



the formula (13). 



One now verifies that the conditional expectations (13) satisfy 

EfoEs =Emin(s,t)- 
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This is true by construction for s < t. The proof for t < s requires a 
tiny computation. One first gets the factors Xg^^ and Yg^^ by applying 
Eg. Applying next Ej we get an extra contribution to X coming from the 
difference between reordering the operators before and after the projection 
by Et which is 

i<j 

One checks that Xs^^ + Xt;s = X^^'' . Similarly, acting next with we get 
an extra contribution to Y which is 

^t;s = Yl^i^^i■''^^^oih■,s) - P'i{U;t)tJ'j{tj;t)]- 
i,j 

One also checks that Y^^^ + Yt-s = y}^\ 

Lastly we prove that the conditional expectations are right neutral by 
measurable multiplication, i.e. E5[6a] = Es[6]a for a € Ag. It is enough 
to check it for ao = e^"^^(*°)e^o^*^*°) and b a product of operators as above. 
If ao multiplies this product of operators on the right, we first have to 
move it to the left by using the commutation relations, then to use the left 
neutrality and apply E.^ and, finally, to move it back to the right using again 
the commutation relations. The difference between the commutators before 
and after having applied E^ is: 

j 

+ 1] (Ao/^ife;.)[^(*o),€^(ts)] +/^oMife>)[e^(^o),4(^.)])• 
J 

This vanishes thanks to the relations l^j[C{to),^'^{tj)] = l^j{tj-s)[C{to), C^ts)], 
and their complex conjugate, valid for to < s < tj. 
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